A  METHOD  FOR  DETERMINING  THE  VIRTUAL  MASS  DISTRIBUTION  AROUND  A  VIBRATING  billF 


KAVAL  SHIP  RESEARCH  AND  DEVELOPMENT  CENTER 


Washington  .D.C.  20007 


A  METHOD  FOR  DETERMINING  THE  VIRTUAL  MASS 
DISTRIBUTION  AROUND  A  VIBRATING  SHIP 


by 

A.N.  Hicks 


This  document  has  been  ‘approved 
for  public  release  and  sale;  its 
distribution  is  unlimited. 


DEPARTMENT  OF  STRUCTURAL  MECHANICS 
RESEARCH  AND  DEVELOPMENT  REPORT 


January  1970 


L 


D  C, 

dT" 

it'  7 

ktsenns 

-  c 


Report  3272 


1 


J 


3n 


The  Nava!  Snip  Research  and  Development  Center  is  a  U.S.  Navy  center  for  laboratory 
effort  directed  at  achieving  improved  sea  and  air  vehicles.  It  was  formed  in  March  1987  by- 
merging  the  David  Taylor  Model  Basin  at  Carderock,  Maryland  and  the  Marine  Engineering 
Laboratory  at  Annapolis,  Maryland.  The  Mine  Defense  Laboratory,  Panama  City,  Florida 
became  part  of  the  Center  in  November  1987. 

Naval  Ship  Research  and  Development  Center 
Washington,  D.C.  20007 


NOW-NSRDC  3960/ 43  ( IP-67) 


DEPARTMENT  OF  THE  NAVY 
NAVAL  SHIP  RESEARCH  AND  DEVELOPMENT  CENTER 

WASHINGTON.  D.  C.  20007 


A  METHOD  FOR  DETERMINING  THE  VIRTUAL  MASS 
DISTRIBUTION  AROUND  A  VIBRATING  SHIP 


by 


A.N.  Hicks 


This  document  has  been  approved 
for  public  release  and  sale;  its 
distribution  is  unlimited. 


January  1970 


Report  3272 


TABLE  OF  CONTENTS 


Page 

ABSTRACT  . 1 

ADMINISTRATIVE  INFORMATION  . 1 

INTRODUCTION  .  1 

NATURE  OF  THE  MODIFICATION  REQUIRED  .  3 

MATHEMATICAL  MODEL  FOR  AN  AXISYFNETRIC  SHIP  . . .  4 

COMPARISON  OF  RESULTS  WITH  KNOWN  EXACT  FORCE  DISTRIBUTIONS  .  8 

VIBRATING  PROLATE  SPHEROID  . 8 

INFINITE  CIRCULAR  CYLINDER  .  10 

APPLICATION  OF  THE  TECHNIQUE  TO  NCNAX I SYMMETRIC  SHIPS  .  12 

APPLICATION  TO  SHIP  VIBRATION  .  13 

COMPARISON  WITH  FULL-SCALE  SHIP  VIBRATION  RESULTS  .  IS 

USE  OF  THE  FULL  BOUNDARY  EQUATION  .  17 

CONCLUSIONS  . 21 

APPENDIX  -  EXACT  SOLUTION  FOR  TRANSVERSE  VIBRATION  OF  AN  ELLIPSOID 

OF  REVOLUTION  .  27 

REFERENCES  .  30 

LIST  OF  FIGURES 

Page 

Figure  1  -  Geometry  and  Coordinate  System  .  23 

Figure  2  -  Coordinate  System  for  Vibrating  Ellipsoid  .  23 

Figure  3  -  Comparison  of  Three-  and  Two-Dimensional  Flow  Approxi¬ 
mations  with  the  Exact  Solution  for  a  Vibrating 
Ellipsoid  (L/B=10)  24 

Figure  4  -  Comparison  of  Three-  and  Two-Dimensional  Flow  Approxi¬ 
mations  with  the  Exact  Solution  for  a  Vibrating 
Ellipsoid  (L/B=5)  25 

Figure  5  -  Measured  and  Computed  Mode  Shapes  for  a  Destroyer  .  26 

LIST  OF  TABLES 

Page 

Table  1  -  Coefficients  for  the  Ellipsoid  Vibration  Shapes  .  9 

Table  2  -  Values  of  f(z)/*pb  vq  for  an  Infinite  Cylinder  .  11 

Table  3  -  Measured  and  Computed  Frequencies  for  a  Destroyer  .  15 

ii 


ABSTRACT 


A  method  is  presented  for  determining  the  three- 
dimensional  virtual  mass  distribution  associated  with  the 
vertical  girder  vibration  of  ships.  The  method  was  de¬ 
veloped  for  use  with  a  lumped  mass/weightless  beam  ship 
representation  and  is  based  on  a  set  of  dipole  distri¬ 
butions  along  the  ship  axis.  It  provides  a  virtual  mass 
matrix  with  off-diagonal  elements  and  enables  all  the  vi¬ 
bration  frequencies  and  shapes  of  the  ship  to  be  computed 
from  a  single  matrix  equation.  The  usual  method  for 
determining  the  frequencies  and  shapes  uses  a  separate 
mass  matrix  for  each  mode.  The  method  is  preferable  to 
the  standard  one  for  short  or  unusual  ships  and  mode 
shapes,  or  where  it  is  desirable  to  include  all  modes  in 
a  single  equation.  However,  if  separate  consideration 
of  each  mode  is  acceptable,  the  standard  technique  is 
simpler  for  normal  ships. 

ADMIN I STRAT I VE  I NFORMAT I ON 


The  author  of  this  report  is  a  Senior  Scientific  Officer  of  the 
Naval  Construction  Research  Establishment  at  Dunfermline,  Fife,  Scotland. 
During  1969,  he  was  assigned  duties  as  exchange  scientist  in  the  Ship  Pro¬ 
tection  Division  of  the  Department  of  Structural  Mechanics,  Naval  Ship 
Research  and  Development  Center;  his  salary  and  expenses  during  this 
assignment  were  paid  by  the  United  Kingdom.  The  work  reported  herein  was 
performed  under  Naval  Ordnance  Systems  Command  Task  UF17-554-504 ,  with 
funding  support  only  for  computer  time  and  printing  of  the  report. 

INTRODUCTION 


In  problems  concerned  with  the  vibration  of  ships,  it  has  long  been 
appreciated  that  the  effect  of  the  surrounding  water  must  be  allowed  for 
if  reasonably  accurate  predictions  are  to  be  made  of  the  frequencies  of 
vibration  and  of  the  mode  shapes.  At  almost  all  frequencies  of  interest, 
the  principal  effect  of  the  water  is  a  very  large  increase  in  the  effective 
inertial  mass  of  the  ship;  at  heaving  but  at  pitching  frequencies,  there  is 
also  a  considerable  damping  effect  due  to  the  generation  of  surface  waves. 
For  many  ships,  the  additional  mass  (virtual  mass)  due  to  the  water  is 
between  1  1/2  and  twice  the  total  displacement  of  the  ship  and  its  effect 
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is  consequently  very  significant.  Todd  gives  an  historical  survey  of 

2 

the  study  of  virtual  mass  and  its  application  to  ship  vibration.  Kaplan 
gives  a  comprehensive  critical  review  of  the  literature  of  the  subject. 

Most  current  methods  of  allowing  for  the  virtual  mass  are  adap¬ 
tations  of  a  technique  introduced  very  early  by  Lewis. *  He  noted  that 
because  of  the  long  slender  nature  of  ships,  the  fluid  flow  around  the 
ship  would  be  largely  confined  between  transverse  planes  because  of  its 
transverse  motion  and  so  could  be  approximated  at  each  ship  cross  section 
by  the  two-dimensional  flow  aiound  an  infinite  cylinder  of  the  same  shape 
as  the  cross  section.  Such  two-dimensional  flows  are  fairly  easy  to 
determine.  The  method  has  since  been  called  strip  theory.  To  allow  for 
the  existence  of  some  flow  parallel  to  the  ship  axis,  Lewis  computed  the 
exact  kinetic  energy  for  vibrating  ellipsoids  of  revolution  and  also  the 
kinetic  energies  for  the  same  ellipsoids  of  revolution  under  the  as¬ 
sumptions  of  strip  theory.  For  any  particular  ship,  the  virtual  mass 
distribution  along  the  ship  deduced  from  the  two-dimensional  flow  so¬ 
lutions  could  then  be  reduced  in  the  ratio  of  the  two  kinetic  energies 
calculated  for  an  ellipsoid  of  the  same  length/beam  ratio.  This  is  still 
the  standard  method  of  computing  the  virtual  mass  distribution  and  the 

two-dimensional  flows  have  been  calculated  for  a  much  greater  variety  of 
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cross-sectional  shapes  than  were  given  by  Lewis.  However,  Tayloi  showed 
that  Lewis  had  not  used  the  best  possible  boundary  condition  for  his 
exact  ellipsoid  solution.  Taylor  presented  a  different  set  of  reduction 
factors  to  allow  for  the  three-dimensional  flow  effect.  He  also  gave  a 
third  set  of  reduction  factors  deduced  from  the  transverse  vibration  of  an 
infinite  circular  cylinder  with  a  sinusoidal  distribution  in  the  vibration 
amplitude  along  its  length.  Most  vibration  calculations  have  used  his 
ellipsoid  correction  factor  but  at  least  one  method^  uses  the  cylinder 
factor. 

The  above  technique  works  very  successfully  for  the  oasic  two  or 
three  vibration  modes,  but  it  has  several  disadvantages.  First,  in 
addition  to  the  variety  of  three-dimensional  flow  correction  factors 
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available,  each  type  of  reduction  factor  depends  on  the  particular  type 
of  motion;  a  different  factor  is  required  for  each  vibration  mode.  Lewis 
gave  correction  factors  for  heaving,  pitching,  and  two-  and  three-node 
vertical  vibration,  Taylor  gave  a  correction  factor  only  for  two-node 
vibration.  Since  most  methods  of  computing  the  vibration  frequencies  in¬ 
volve  eigenvalues  of  either  a  differential  equation  or  a  matrix,  this 
means  that  a  different  equation  or  matrix  must  be  solved  for  each  mode; 
this  requirement  considerably  complicates  programming  a  computer  to  carry 
out  the  task.  Moreover,  for  some  shapes  (as  discussed  at  some  length  in 
Reference  2),  there  is  the  possibility  that  the  correction  factor  may  vary 
along  the  length  of  the  ship,  the  reduction  possibly  being  greater  near 
the  ends  than  at  the  center.  There  is  no  present  means  of  allowing  for 
such  a  variation.  Finally,  in  one  particular  type  of  problem,  namely  the 
hull  whipping  induced  by  underwater  explosions,  in  order  to  be  able  to 
deal  with  nonlinear  effects  conveniently  it  is  necessary  to  represent  the 
equations  of  motion  of  the  ship  as  a  single  matrix  equation  including  the 
effects  of  all  modes  simultaneously.  The  present  type  of  correction 
factor,  varying  from  mode  to  mode,  excludes  such  a  representation.  Even 
in  the  completely  linear  case,  the  vibrations  of  a  damaged  ship  involve 
mode  shapes  for  which  no  correction  factors  are  available.  These  dis¬ 
advantages  prompted  the  presen.  attempt  to  find  an  alternative  method  of 
allowing  for  the  virtual  mass  effect. 

NATURE  OF  THE  MODIFICATION  REQUIRED 

The  present  "strip  theory"  described  above  relates  the  force  on 

each  transverse  section  of  the  ship  to  the  motion  of  that  section  alone, 

the  overall  (three-dimensional)  correction  factor  being  constant  along 

the  length.  If  the  ship  is  divided  into  n  sections,  with  n  large  enough 

that  the  variation  in  cross-sectional  shape  is  small  in  each  section, 

t"h 

then  the  vertical  hydrodynamic  force  F^  on  the  i  '*  section  is  related  to 
the  average  vertical  displacement  v.  of  the  section  by 

F  .  =  -  m  .  v. 
wi  wi  ' 1 

where  dots  denote  differentiation  with  respect  to  time  and  m  .  is  the 

wi 

virtual  mass  of  the  section.  The  equation  for  the  forces  on  all  sections 


of  the  ship  may  then  be  written  in  the  matrix  form 


where  F  and  y  are  column  vectors  and  M  is  a  diagonal  matrix.  This 
_w  i  w 

standard  strip  theory  equation  clearly  cannot  be  valid  for  general  motions. 
If  one  section  of  the  ship  is  accelerating  upward,  it  will  certainly  ex¬ 
perience  a  downward  hydrodynamic  force  and  so  diagonal  elements  in  M  are 
necessary.  However  it  will  also  induce  a  downward  fluid  acceleration 
around  it  as  fluid  moves  to  fill  the  space  it  is  vacating.  This  flow  will 
be  around  neighboring  sections  and  so  will  induce  a  downward  force  on  these 
too.  This  means  that  for  Equation  [1]  to  be  true,  M  should  contain  off- 
diagonal  elements.  In  fact,  all  elements  of  M  will  be  nonzero,  but  since 
each  section  will  principally  affect  its  nearest  neighbors,  the  magnitude 
of  the  elements  will  decay  rapidly  away  from  the  main  diagonal.  The 
problem  is  to  determine  the  elements  of  M  .  Since  the  general  ship  prob¬ 
lem  is  very  difficult,  attention  is  directed  first  to  the  case  of  an 
axisymmetric  ship. 


MATHEMATICAL  MODEL  FOR  AN  AXISYMMETRIC  SHIP 


For  such  a  ship,  it  is  possible  to  satisfy  fairly  well  the  fluid 
boundary  condition  of  equality  of  hull  and  fluid  velocities  along  the 
normals  to  the  hull  by  means  of  a  distribution  of  vertical  dipoles  dis¬ 
tributed  along  the  axis  of  the  ship.  This  distribution  should  normally 
be  continuous,  but  if  the  ship  is  considered  to  be  divided  into  a  number 
of  sections  and  the  number  is  large  enough,  then  it  should  be  reasonable 
to  assume  that  the  line  distribution  in  each  section  i  has  a  constant 
strength  ik. 

Since  the  velocity  potential  at  the  point  (r,9,z)  due  to  a  dipole 
at  the  point  (o,o,s)  (referred  to  the  cylindrical  polar  coordinates  of 
Figure  1)  is  given  by 


where  u  =  -  V4>  , 
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the  potential  due  to  a  line  distribution  of  dipoles  of  strength  u.  per 
unit  length  and  extending  from  z ^  ^  t0  zi+i/2 

/  i+1/2 
I  ds 

6.  (x,0,z)  =  u-  r  cos  0 - xn 

J  r  2  ,  ^2,  7 

zi-l/2  lr  *(2-s)  1 


whence 


ik  cos  9 


z-z .  , 
1-1/2 


'i+1/2 


,  7  1/2  ,  ,  1/2 
[r  +(?_2i-l/2)  1  [r+(2-z.+1/2)2] 


'i+1/2 


The  radial  velocity  u^  produced  by  the  distribution  is  given  by 


3*.  u1  cos  e  j  z-zi+l/2 

Uri'Z  3r  ~  2  |  .  1/2  "  . 

|tr  »(2-ii!/2)]  [r  *(2-zltl/2)  ] 


r  (z-zi-l/2> 


r  (z-zitl/2) 


[r2-Cz-ii_l/2}2l  [r2«(i-3,,1/2)2) 


and  the  longitudinal  velocity  uz  by 


* ;  (  i 

uiiCz)  ’-inr=  “i r  cos  8  {— — 

Ur  *(i- 


.  2 

zi*l/2>  I  lr 


*‘Z-Zi-1/2S  1 


3/2  l"! 


The  total  velocities  produced  by  all  the  distributions  are  therefore 


UrW  =Z)  Urj(Z)  a"d  U2Cz)  “I^u2j(z) 


The  boundary  condition  at  the  hull  is  that  the  velocities  of  the 
fluid  and  the  hull,  resolved  in  the  direction  of  the  normal  to  the  hull. 
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should  be  equal .  Ft  the  axisymmetric  shape  considered,  the  angle  X  be¬ 
tween  the  normal  to  the  hull  and  the  radial  direction  is  given  by 
X  ~  -  tan-*  b'  where  b(z)  is  the  radius  of  the  ship  and  b'  =  db/dz.  The 
boundary  condition  is,  therefore, 

v  cos  X  -  v'  b  sin  X  =  u  cos  X  +  u  sin  X  f£i 

r  z  1  1 

where  v(z)  is  the  distribution  of  vertical  velocity  along  the  length.  This 
condition  allows  for  both  shearing  and  flexing  of  the  ship.  For  the  more 
interesting,  lower  frequency  modes  of  vibration,  v'b  will  be  small  com¬ 
pared  to  v  since  the  wavelength  will  be  much  greater  than  the  half  beam. 

X  will  also  normally  be  small  except  possibly  in  the  immediate  neighbor¬ 
hood  of  the  stern.  The  term  V'b  sin  X  will  therefore  be  very  small.  The 
flow  along  the  ship,  (up  is  produced,  partly  by  the  variation  in  v  along 
the  ship  (i.e.,  by  v')  and  partly  by  the  changing  shape  of  the  ship  (b') . 

It  too  will,  therefore,  normally  be  small  when  v'  and  b'  are  small.  It  is 
generally  possible  then  to  simplify  the  boundary  condition  (Equation  [5]) 
to 

v  =  ur  [6] 

except  when  rapid  changes  occur  in  either  the  ship  underwater  cross- 
sectional  area  or  in  the  velocity  distribution  along  its  length  (e.g.,  for 
high  modes).  This  simplification  in  the  boundary  condition  slightly 
reduces  the  amount  of  data  needed  to  specify  the  ship  and  its  motion.  It 
is  roughly  equivalent  to  neglecting  rotary  inertia  in  the  dynamic  equations 
of  the  ship  itself. 

Neither  of  the  boundary  conditions,  Equations  [5]  and  [6],  can  be 
satisfied  everywhere  by  the  assumed  velocity  potential.  However,  either 
can  be  satisfied  at  up  to  n  "collocation”  points  along  the  ship  and  the 
boundary  condition  will  not  be  seriously  violated  anywhere  if  these  points 
are  suitably  chosen.  The  most  convenient  choice  for  the  collocation 
points  is  at  the  midpoints  of  the  ship  sections,  and  satisfying 
Equation  [6]  at  these  points  gives 

y  =  -2  A  y  [7] 

Here  b  is  the  maximum  value  of  b(z),  y  (z)  is  the  vertical  velocity  of  the 
ship,  and  the  matrix  A  =  (a.^)  where 
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a.  .  = 
13 


(Zi"2j+V2) 


b. 


,2 


1/2 


i  \  tbi  +(zi-zj'+£j/2)  ] 


( z.-Zj-y2 ) 
[bAcz.-z.-*^}2] 


1/2 


bi  CZi-Zj4£j/2) 


b.  (z. -z .-9.-/2) 
1  1  J  J 


2  2  ^2  2  2  ^2 
[b.S(zi-z.+il./2)Z]  [b.S(z.-z .  -l.nt) 


IS] 


Equation  [7]  may  then  be  inverted  to  determine  the  strengths  p  of  the  line 
distributions  necessary  to  satisfy  the  boundary  condition.  Equation  [6] , 
for  the  given  velocity  distribution  y.  This  gives 


,2  A"1  • 

p  =  b  A  y 


The  upward  vertical  force  per  unit  length  on  the  section  at  the 
point  z,  is  given  by 

7T 

f.  =  -  2  [  p  cos  0  b.  d6 

1  V 

where  p  =  p4>  is  the  fluid  pressure.  Thus 
f..  =-  2p  J*  6  cos  6  d6  =-  2p^ 


i  n 
b 
l 


;2 

j  =  l  o 


cos  9  d0j 


n 


£-J  13  3 

1  3  =  1 


where 


eii  ’  b- 


z.-z.+£./2)‘ 
13  3 


z.-z. -?../2) 
13  3 


2  7  1/2  “7  7 

U*>i  +CZ.-Z.+4./2)2]  [b.  +(z.-z.-il./2)/] 


1/2 


[9] 


Since  in  any  actual  motion,  conditions  along  the  ship  will  be  continuous, 
the  upward  force/unit  length  at  the  center  of  the  section  will  approximate 
the  average  force/unit  length  over  the  whole  section  and  the  total  upward 
force  on  the  section  will  be  approximately 


7 


that  is 


[10] 


F. 

l 


£.  f. 


X  1 


n 


Uj 


•  3  _1  *** 

F  =  -  upbB  )i  =  -  ub  p  B  A  y 

where  B  is  the  matrix  (3..).  The  required  inertial  water  mass  matrix  M 

1 3  n  w 

is  therefore  given  by 


M„  =  7tb3p  B  A-1  [111 

w  1  i 

This  matrix  depends  only  on  the  shape  of  the  ship  and  is  completely  inde¬ 
pendent  of  the  type  of  motion  (or  vibration  mode) .  With  M  known,  the 
force  distribution  for  any  vertical  acceleration  distribution  y  is  readily 
found.  A  short  computer  routine  has  been  written  to  compute  the  non- 
dimensional  matrix  BA-1  and  the  force  distributions  deduced  from  it,  by 
Equation  [10],  for  given  distributions  of  vertical  acceleration. 


COMPARISON  OF  RESULTS  WITH  KNOWN  EXACT  FORCE  DISTRIBUTIONS 


Exact  solutions  are  known  for  two  forms  of  fluid  flow  that  are 
suitable  for  comparison,  namely,  the  flows  around  a  vibrating  prolate 
spheroid  and  around  an  infinite  circular  cylinder  whose  transverse 
velocity  varies  sinusoidally  along  its  length.  The  exact  results  for 
both  cases  have  been  compared  with  the  results  from  the  foregoing  analysis. 

VIBRATING  PROLATE  SPHEROID 

The  methods  involved  in  the  solution  of  the  flow  about  a  prolate 
spheroid  are  discussed  in  some  detail  by  Lamb,  (see  page  139  of  Reference 

3 

6) .  Lewis  gave  the  first  solution  in  connection  with  transverse  shear 

4 

vibrations  and  Taylor  gave  a  second  solution  using  a  different,  more 

realistic  boundary  condition  involving  both  flexure  and  shear.  This  type 

7 

of  motion  has  also  since  been  investigated  by  Landweber  and  Macagno. 

In  principle,  the  vibrating  ellipsoid  can  be  solved  exactly  for 
any  arbitrary  transverse  velocity  distribution  but  in  practice  only  dis¬ 
tributions  represented  by  low  order  polynomials  are  required.  Using  the 
analysis  outlined  in  the  Appendix,  a  short  computer  routine  was  written 
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to  evaluate  the  force  distribution  on  an  ellipsoid  with  a  velocity  distri¬ 
bution  representable  by 

S 

v(2)  -JT vnCz/a)n-' 
n=l 

where  z  is  the  distance  along  the  axis  of  symmetry  from  the  center  of  the 
ellipsoid  of  length  2a  (see  Figure  2).  This  velocity  distribution  is 
sufficient  to  approximate  heaving,  pitching,  and  the  first  three  whipping 
modes  of  ship  vibration.  The  boundary  condition  used  in  the  solution  is 
that  of  Taylor  and  allows  for  flexure  as  well  as  shear. 

Figures  3  and  4  compare  the  transverse  force  distributions  given 
by  the  approximate  analysis  with  the  values  given  by  the  exact  solution 
and  also  with  the  values  given  by  the  strip  method,  using  the  Lewis 
correction  factors  since  these  are  available  for  four  of  the  five  modes. 
The  values  used  for  the  coefficients  (v^,...,v^)  for  these  cases  are 
given  in  Table  1. 

TABLE  1 


Coefficients  for  the  Ellipsoid  / ibration  Shapes 


^"'^Coefficient 

Mode 

V1 

V2 

V, 

V4 

vs 

Heave 

1 

0 

0 

0 

0 

Pitch 

0 

1 

0 

0 

0 

2-node  Vertical 

-0.200 

0 

1 

0 

0 

3-node  Vertical 

0 

-0.429 

0 

1 

0 

4-node  Vertical 

0.0274 

0 

-0.534 

0 

1 

The  first  four  mode  shapes  are  those  used  by  Lewis,  although  the  exact 
analysis  used  the  better  Taylor  boundary  condition.  The  Lewis  correction 
factor  and  the  approximate  three-dimensional  analysis  both  use  the  shear 
type  boundary  condition,  and  the  three-dimensional  analysis  also  assumes 
that  the  rate  of  change  of  the  radius  along  the  length  is  small.  The 
fifth  mode  shape  has  been  chosen  to  have  nodes  at  0.155L,  0.38L,  0.62L, 
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and  0.845L,  where  L  is  the  total  length.  These  are  about  the  correct 
positions  for  destroyers,  but  the  resulting  shape  gives  rather  too  much 
prominence  to  the  ends  and  too  little  to  the  central  section. 

For  the  L/B=10  ellipsoid,  the  results  for  both  the  three-dimensional 
and  the  two-dimensional  (strip  theory)  approximations  agree  well  with  the 
exact  analysis  and  there  is  little  to  choose  between  them.  The  three- 
dimensional  approximation  is  slightly  better  near  the  center  of  the 
ellipsoid  where  changes  in  the  radius  are  smallest,  but  the  two-dimensional 
approximation  is  better  at  the  ends  where  the  radius  is  changing  rapidly. 

The  results  for  the  L/B=5  ellipsoid  are  very  similar  but  the  di¬ 
vergence  from  the  exact  solution  is  quite  serious  near  the  ends  for  both 
approximations  for  modes  as  low  as  the  second  vibration  mode.  Once  again, 
there  is  very  little  to  choose  between  the  two  approximate  methods. 

For  most  ships  in  which  vibration  frequencies  are  particularly  im¬ 
portant,  the  L/B  ratio  is  near  10  and  both  the  strip  theory  and  the  new 
three-dimensional  approximation  should  give  good  results.  The  divergence 
near  the  ends  for  the  lower  L/B  ratio  is  caused  by  using  the  approximate 
form  of  the  boundary  condition.  Equation  [6].  From  this  point  of  view, 
the  ellipsoid  is  rather  a  poor  shape  since  the  radius  changes  extremely 
rapidly  near  the  ends.  The  radius  changes  are  much  less  severe  for 
typical  ship  forms,  and  either  approximate  method  would  give  better 
results.  Use  of  the  exact  boundary  condition.  Equation  [5],  would  be 
simple  in  the  case  of  the  ellipsoid  because  the  velocity  distribution  and 
rate  of  change  of  radius  are  easily  defined,  but  is  scarcely  worth  the 
effort  for  ship  forms  (unless  rotary  inertias  are  being  considered) .  This 
point  is  considered  later. 

INFINITE  CIRCULAR  CYLINDER 

4 

Taylor  was  the  first  to  consider  this  case.  The  infinitely  long 
circular  cylinder  was  assumed  to  have  a  transverse  velocity  distribution 

v(z)  =  vq  cos  kz 

Taylor  gave  the  velocity  potential  for  this  distribution  as 


cos  9  cos  kz 


Vo 

$  = - = - 

k  K'(kb) 

where  b  is  the  cylinder  radius  and  is  a  modified  Bessel  function  of 
the  second  kind.  The  force  distribution  may  be  found  as  before  and  is 

f  f  .  .7  V  Vkb) 

f(z)  =-2  I  p  cos  6  b  d9  =-  2pb  |  cos  0  d8  =-  apb“ - - - cos  kz 

J  *  J  kb  K.'(kb) 

o  o  lv  J 

(force/unit  length) 

The  wavelength  X  of  the  velocity  distribution  is  X  =  2a/k, 

Since  the  approximate  analysis  is  based  on  a  body  of  finite  length, 
it  cannot  give  a  uniformly  good  representation  of  the  infinite  cylinder. 
However,  if  it  is  used  to  represent  three  complete  wavelengths  of  the 
cylinder,  the  flow  in  the  central  wavelength  should  be  approximately 
correct.  Kith  three  wavelengths,  the  program  restricts  the  number  of 
sections  in  each  wavelength  to  six.  Table  2  compares  the  results  from  the 
approximate  analysis  with  the  exact  results.  In  this  case,  the  cylinder 
is  of  uniform  diameter  so  that  the  boundary  conditions,  Equations  [5]  and 
[6],  in  the  three-dimensional  approximation  are  equivalent.  Inaccuracies 
in  the  solution  are  due  either  to  the  coarseness  of  the  representation  or 
to  the  finite  length  of  the  cylinder  in  che  three-dimen:  lonal  approxi¬ 
mation. 


TABLE  2 

2  # 

Values  of  f(z)/Ttpb  vq  for  an  Infinite  Cylinder 


^X/b 

a/6 

~/2 

m 

3r/2 

11  -/6 

13  -/6 

5-/2 

17  -/6 

5 

exact 

-0.448 

0 

jfMfH 

0 

-0.448 

9 

0.448 

3-0 

-0.0002 

B§SI119 

mm 

-0.0014 

0.383 

10 

exact 

■ 

0 

(SB 

0.631 

0 

-0.631 

-0.631 

0 

mm 

3-D 

H 

-0.0003 

B 

0.609 

-0.0013 

-0.611 

-0.616 

-0.0189 

HI 

15 

exact 

-0.762 

0 

IB 

IB 

0 

-0.762 

-0.762 

0 

0.762 

3-0 

-0.741 

-0.0001 

B 

B 

-0.0004 

-0.744 

-0.745 

-0.0096 

0.676 

20 

exact 

SB 

0 

B 

0.826 

0 

-0.826 

-0.826 

0 

0.826 

3-D 

EH 

-0.0007 

B 

0.339 

0.0007 

-0.839 

-0.840 

-0.0015 

0.823 

For  all  values  of  A/b,  the  forces  are  cJ early  very  accurate  near  the 
center,  being  within  4  percent  of  the  exact  value  in  all  cases.  The 
points  given  were  the  only  ones  used,  indicating  the  coarseness  of  the 
representation;  it  is  equivalent  in  a  ship  representation  to  using  20 
points  to  represent  the  fifth  vibration  mode  (six  nodes) .  Since  the 
values  from  the  three-dimensional  approximation  are  constant  to  at  least 
z/A  =  7«/6,  the  differences  from  the  exact  result  are  attributable  to  the 
coarse  mesh  rather  than  to  the  finite  length.  For  the  A/b  *s  5  case,  the 
lengths  of  each  section  are  nearly  equal  to  their  radii  whereas  for  the 
A/b  =  20  case,  the  sections  are  nearly  nine  radii  long. 

APPLICATION  OF  THE  TECHNIQUE  TO  NON AXI SYMMETRI C  SHIPS 

Clearly  it  would  be  possible  in  principle  to  extend  the  three- 
dimensional  approximation  by  adding  distributions  of  higher  multipoles 
along  the  ship  axis  and  determining  their  strength  by  satisfying  a 
boundary  condition  around  the  circumference  of  the  ship  as  well  as  along 
its  axis.  This  would,  however,  require  very  large  amounts  of  data  to 
represent  the  ship  shape  as  well  as  the  inversion  of  a  very  large  matrix. 
The  success  of  the  strip  method  using  Lewis  sections  for  the  cross- 
sectional  shapes  points  to  simpler  approaches. 

In  the  solution  of  the  two-dimensional  flows  about  ship-type  cross 
sections,  although  the  velocity  potential  may  in  fact  consist  of  a  super¬ 
position  of  two-dimensional  multipoles  of  all  orders,  the  added  mass  of 
each  section  depends  only  on  the  dipole  term,  the  shape  of  the  section 
determining  its  strength.  In  the  three-dimensional  case,  it  should, 
therefore,  be  approximately  correct  to  account  for  the  shape  of  the 
section  via  the  strengths  of  the  dipole  distributions. 

Since  the  original  work  by  Lewis,  it  has  been  customary  to  represent 
the  added  mass  per  unit  length  of  the  two-dimensional  cross  sections  in 
the  form 


71  >2  r 

7P  b  C 

where  b  is  the  half  beam  of  the  section  and  C  is  a  constant  depending  on 
the  section  shape.  Values  for  C  have  been  computed  for  a  great  variety  of 
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shapes. ^*4,8-12  ^he  above  value  for  the  added  mass  is  also  that  due  to  a 
circular  cylinder  of  radius 

b  .  =  bl/C"  [121 

For  a  r.onaxi  symmetric  ship  then,  each  cross  section  may  be  compared  with 
the  known  shapes  and  C  determined.  Then  Equation  [12]  gives  the  appro¬ 
priate  radius  for  an  axisymmetric  approximation  to  the  actual  shape.  There 
are  no  reasonably  simple  nonaxi symmetric  three-dimensional  flows  with 
exact  solutions  which  can  be  used  for  comparison,  but  the  procedure  should 
give  reasonable  results.  Certainly  in  the  2-node  vibration  mode,  in  the 
central  section  of  the  ship  where  the  added  mass  is  most  important,  the 
technique  will  give  very  good  answers  since  the  flow  in  this  region  is 
very  nearly  two-dimensional  and  the  method  is  exact  in  the  two-dimensional 
case. 


APPLICATION  TO  SHIP  VIBRATION 
5  IS 

Most  current  techniques  *  for  the  determination  of  ship  natural 
frequencies  by  purely  theoretical  means  depend  on  finite-element  lumped 
mass  approaches.  These  represent  the  ship  as  a  series  of  lumped  masses 
interconnected  by  weightless  elastic  beams.  All  applied  forces,  including 
distributed  inertial  forces  and  moments,  are  approximated  by  equivalent 
point  forces  and  moments  applied  to  the  lumped  masses.  It  is  then 
possible  to  compute  a  stiffness  matrix  K  such  that  when  no  moments  are 
applied,  the  forces  F  required  at  the  masses  to  statically  maintain  a 
given  displacement  shape  y  are  given  by 

F  =  K  y 

Neglecting  buoyancy  forces  and  rotary  inertia  (these  can  easily  be  included 
if  desired),  the  only  forces  on  a  ship  in  still  water  are  inertial  forces 
Fj  and  hydrodynamic  forces  F^.  If  the  values  of  the  lumped  masses  are 
- • • > 

Fj  =  -  M  y 
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where  M  is  the  diagonal  matrix  with  elements  (ra^) .  The  hydrodynamic 
forces  are  given  by  Equations  [10],  i.e., 

F ,  =  -  M  y 
_w  w  i. 

with  as  given  in  Equation  [11]. 

The  equation  of  ship  motion  is  then 

K  y  =  F  =  -  tr  -  M.  y  [13] 

i.e.  (M  +  M  )  y  +  K  y  =  0 

v  -  _ 

The  natural  vibration  frequencies  are  the  eigenvalues  of  this  matrix 

equation  and  the  mode  shapes  are  the  corresponding  vectors.  In  the  strip 

method  for  the  hydrodynamic  flow,  M  and  M  are  both  diagonal  matrices. 

w 

Since  K  is  symmetric,  the  equation  can  then  be  transformed  into 

z  +  S  z  =  0 

where  S  =  (M  +  M  )-1/2  K  (M  +  M  )“1/2  and  z  =  (M  +  M  )1/2  y. 

w  w  w  .. 

S  is  symmetric  and  its  eigenvalues  are  easily  found  by  any  of  the  standard 
routines  for  eigenvalues  and  vectors  of  symmetric  matrixes.  For  the  strip 
method,  however,  Mw  depends  on  the  mode  shape  being  investigated  and  a 
different  matrix  S  must  be  used  for  each  mode. 

In  the  proposed  three-dimensional  analysis,  the  matrix  Mw  is  found 
as  a  full  matrix  with  a  dominant  diagonal  but  no  zero  elements.  Equation 
[13]  may  then  be  written 

y  +  S.  y  =  0  ;  S.  »  (H  +  M  )‘X  K 

_  1  -  I  iV 

and  the  eigenvalues  found  directly.  will  not,  however,  be  symmetric, 
and  this  restricts  the  available  range  of  computer  routines.  The  full  M 
matrix  is  actually  not  symmetric  but  the  degree  of  asymmetry  is  not 
large  except  for  extreme  shapes.  It  may  be  nrtifically  made  symmetric  by 
replacing  all  elements  in^  by  1/2  (in^  +  m^..^).  This  procedure  was 
carried  out  for  the  examples  used  to  check  the  three-dimensional  theory  and 
in  no  case  did  it  change  the  resulting  forces  by  more  than  3  percent.  Since 
the  vibration  frequencies  depend,  approximately,  on  the  square  root  of  the 
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mass,  this  difference  is  negligible.  A  standard  Cheieski  decomposition 
14 

routine  may  then  be  used  to  generate  a  matrix  L  such  that 

LLT  =  0?  +  M  ) 

Equation  [13]  then  transforms  into 

z  *■  S2  z  =  0 

-1  -T 

where  S£  =  L  K  L  .  S9  is  now  symmetric  and  its  eigenvalues  and  vectors 
may  again  be  found  by  standard  symmetric  routines. 

COMPARISON  KITH  FULL-SCALE  SHIP  VIBRATION  RESULTS 

Frequencies  and  mode  shapes  for  overall  hull  vibrations  were  measured 
recently  or.  a  World  War  II  2500-ton  destroyer,  HMS  ROEBUCK.  Table  3  com¬ 
pares  measured  frequencies  with  those  calculated  using  both  standard  strip 
theory  and  the  three-dimensional  flow  approximation.  The  strip  theory 
results  were  obtained  by  using  Lewis  three-dimensional  correction  factors 
since  these  are  available  for  heave,  pitch,  and  thp  first  two  vibration 
modes  and  could  be  estimated  for  the  third  and  fourth  modes  by  extrapo¬ 
lation.  Also  included  in  the  table  are  the  results  calculated  for  the 
three-dimensional  flow  approximation  using  the  full  boundary  condition. 
Equation  [5],  as  described  later. 


TABLE  3 

Measured  and  Computed  Frequencies  for  a  Destroyer 
(Frequencies  are  given  in  hertz) 


Computed 

Mode 

Measured 

Strip  Method 

3-D  Flow  Using 
Equation  [6] 

3-D  Flow  Using 
Equation  [5] 

-- 

0.20 

0.22 

0.20 

0.25 

0.20 

0.25 

First  Vib 

1.68 

)  ■  o4 

1.64 

1.66 

Second  Vib 

3.35 

3.16 

3.20 

3.22 

Third  Vib 

4.97 

4.92 

4.96 

5.00 

Fourth  Vib 

6.63 

6.90 

6.96 

7.03 

The  results  of  all  three  methods  are  clearly  fairly  reasonable  and 
there  is  veiy  little  to  choose  between  them.  The  three-dimensional  flow 
theory  using  the  approximate  boundary  condition*  Equation  [6] ,  gives 
results  very  similar  to  the  strip  method  using  Lewis  correction  factors. 
Since  the  approximate  boundary  condition.  Equation  [6] ,  is  equivalent  to 
the  Lewis  original  one,  the  agreement  really  is  a  justification  of  the  use 
of  strip  theory.  The  results  for  the  three-dimensional  flow  with  the  more 
exact  boundary  condition.  Equation  [5] ,  were  slightly  better  and  reduced 
the  error  in  the  first  mode  frequency  from  2  1/2  to  1  1/4  percent.  The 
strip  theory  gave  the  same,  improved,  result  for  the  first  mode  frequency 
when  the  Taylor  rather  than  the  Lewis  correction  factor  was  used.  How¬ 
ever,  Taylor  does  not  give  values  for  the  reduction  factor  for  the  other 
modes.  For  the  strip  theory  results,  a  separate  three-dimensional  factor 
had  to  be  applied  for  each  mr-de. 

The  differences  between  the  computed  and  measured  frequencies,  how¬ 
ever,  were  slightly  larger  than  the  differences  between  the  computed 
values  themselves,  indicating  that  the  remaining  errors  were  probably  due 
to  factors  otner  than  the  hydrodynamics .  There  is  still  some  doubt  over 
the  ship  mass  distribution,  the  material  to  be  included  in  calculating 
the  stiffness  distribution  along  the  ship,  and,  probably  most  important, 
the  best  method  for  calculating  the  shear  area  distribution  along  the  ship. 

The  calculated  mode  shapes  were  practically  identical  for  the  three 
methods;  they  are  compared  with  the  experimental  shapes  for  the  first 
three  vibration  modes  in  Figure  5.  The  difference  between  the  measured 
and  calculated  shapes  for  the  first  two  modes  was  less  than  the  scatter  in 
the  experimental  values,  but  there  was  a  definite  difference  between  the 
shapes  near  the  bow  for  the  third  mode. 

The  computer  program  which  produced  the  results  for  both  the  three- 
dimensional  flow  approximations  is  a  modification  of  part  of  the  FORTRAN 
IV  program  described  in  Reference  13  and  is  run  on  an  IBM  7090  computer. 
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USE  OF  THE  FULL  BOUNDARY  EQUATION 


Where  fairly  sharp  changes  occur  in  the  cross-sectional  shape,  the 
three-dimensional  flow  approximation  can  be  improved  by  using  the  full 
boundary  condition  (Equation  [5])  in  place  of  the  approximate  condition 
(Equation  [6]).  Using  the  finite  element/ lumped  mass  approach,  with  n 
masses,  the  elastic  nature  of  the  ship  may  be  represented  by  the  equation 


F~ 

b” 

m 

.9. 

1/ 

C 

U 

where  ?  and  Q  are  n-component  vectors  that  respectively  represent  the 
forces  and  moments  which  must  be  applied  to  the  lumped  masses  to  maintain 
statically  the  displacements  y  and  bending  rotations  Y  (Y  excludes  shear 
deformations).  A,  B,  and  C  will  be  (n  x  n)  matrixes  whose  elements  are 
given  m  Reference  13.  In  vibration  applications,  the  forces  and  moments 
are  just  the  inertial  and  hydrodynamic  forces  and  moments  acting  on  the 
regions  around  the  lumped  masses.  F  and  Q  will,  therefore,  be  given  by 


■«’-  — 

F 

1 

II 

M 

— 

0 

y 

M  , 
wl 

M  0 
w2 

- ! 

:  : 

_ 1 

_Q_ 

_0 

R_ 

_Y_ 

M  - 

WJ 

M  . 
w4~* 

LIJ 

where  M  and  R  are  diagonal  matrixes  whose  elements  are  the  ship  masses  and 

rotary  inertias  at  the  lumped-mass  oositions.  The  matrixes  M  ,,  M  M  ,, 

r  wl  w2  wo 

and  M  ^  will  be  hydrodynamic  mass  and  inertia  matrixes. 

The  above  elastic  matrices  are  of  order  (2n  x  2n) ,  and  there  are 
several  alternatives  for  the  hydrodynamic  formulation.  The  number  of 
sections  on  which  the  dipole  distributions  are  defined  may  be  doubled  or 
the  same  number  of  distributions  may  be  kept,  but  the  distributions  them¬ 
selves  must  be  given  a  linear  variation  instead  of  being  constant.  With 
each  technique,  either  the  full  or  the  approximate  boundary  conditions 
could  be  used  although  the  full  condition  is  the  more  appropriate.  To 
determine  the  unknown  strengths  of  the  2n  dipole  distributions,  2n  col¬ 
location  points  would  be  required.  As  an  alternative,  the  same  dipole 
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distribution  can  be  assumed  but  be  made  to  satisfy  the  full  boundary  con¬ 
dition.  It  can  be  used  to  calculate  both  the  forces  and  the  moments  at 
the  lumped  masses.  This  latter  technique  is  in  accord  with  using  rotary- 
inertia  terms  in  the  elastic  representation  of  the  ship  instead  of  dou¬ 
bling  the  number  of  elastic  sections  and  lumped  masses,  and  is  adopted 
here. 

th  • 

The  vertical  velocity  at  the  i  mass  is  y^  and  the  angular 
velocity  is  Therefore,  substituting  for  the  angle  A,  the  full  boundary 
condition  [5]  at  the  position  z^  is 


n  n 

y.  +  b.’o.  '  Y.  -  u  .  (z.)  -  b'.  u  .  (z.)  =  -^7  5”* a.  .  11. 

111  rJ  1  1  z r  r  h2^  -1 

j=l  L  J  j=l 


[14] 


where 


a.  .  =  a.  .  +  b  b 
ij  ij 


iv|‘i2,(vyV!)!]  '  -[bi2*(vy 

-V2)J  ) 


Equation  [14]  may  be  written  in  matrix  form  as 


y  +  BY  =  4  A  p  • 

~  b  ~ 

where  B  is  the  diagonal  matrix  with  elements  b^b^'*.  The  strengths  of  the 
dipole  distributions  are  then 

V  =  b2  A'1  y  +  b2  A-1  B  Y  [15] 

Given  the  strengths  of  the  dipole  distributions,  from  Equation  [15],  the 
velocity  potential  is  given  by 

3  =  1 
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where  ^  is  defined  in  Equation  [2] .  This  leads  to  a  vertical  force 
f(zjdz  on  a  length  dz  of  ship;  f(z)  is  given  by 


» -  -2V, ' '  ■  f  ‘"-/1  iJ‘. 

>■'  It  •  >•- w  J  [.’ .  „.r I 

The  hydrodynamic  force  and  moment  on  the  i^1  section  are  therefore 


z  -  z. 


[16] 


vj 


Ji+1/2 


f(z)  dz  and  Q 


ri+1/2 


(z-zp  f(z)  dz 


i-1/2 


Zi-l/2 


whence 


=-  up  /  c. .  p.  ,  Q.  =-  rip  y  d.  .  jj . 

i  jL-j  i]  ]  xi  jLj  ij  3 

j=l  j=l 


In  these,  c„  and  d^  are  given  by 


c.  Jl. 
1 3  i 


■  M'1/2  L  -r,  -r2  „2  )  ♦ 

\  \U  1-  2.  2-j  ^ 


_  ,  /V1/2*rn. 

1/2  l08  U1-i/2<-r1_ 


X2-“2  ,  /V1/2*r2<- 

7~ los  U-l/2«2. 


K) 


Z.-Z  +1./2  z  -z.-9.J2  b.b r/Jl.+x.  b./4.-x,br 

v.h«r-  V  -  _Ji_  J _ a _  „  -  1  3  1  .  K  _  1  1  1  1  .  r  ..  1  l  1  1 

rG  1  *.  ’  X2  i.  ’  al  "  ...2  'B1-“T1 

i  l  1+b^  1+br 

r  2  2 r  2  2~i^2 

r1+  =  (aj+1/2)  +  Bj  ;  r^  =  (^-1/2)*  +  ,  with 

similar  expressions  for  62>  r2+  and  r2  » 

dij/£i2  =  2*  (1+bi2)  [(al+1/2)rl+  "  C«1~l/2)r1_  -  (a2+l/2)r2+  +  (a2-l/2)r2 
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With  these  values  and  uniform  length  sections. 


C  A_1B 


D  A-1B 


[17] 


where  C  and  D  are  the  matrixes  (c.^)  and  (d^),  which  gives  the  hydro- 
dynamic  matrixes  H  j,  Mw2,  M  The  equation  of  ship  vibration  is 
then 


M  +  M  ,  M  - 
wl  w2 

M  -  R  +  M  , 
w3  w4 


=  0 


[18] 


Since  the  rotary  inertia  corrections  are  small,  the  matrixes  M  R,  and 
M  ^  will  be  small.  If  they  are  neglected,  the  equation  for  Y  gives 


Y 

and  the  y  equation  becomes 


y 


T  0-1  n' 

r]  ••  f 

-i 

f] 

LM+Mwl-Mw2C  B 

\y  +  l 

A  -  B  C  B 

J  r-° 

[19] 
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This  equation  is  similar  to  Equation  [14]  but  allows  fully  for  the  effect 
of  changing  cross  sections  and  bending  deformation  in  the  boundary  con¬ 
dition. 

The  force  distributions  given  by  Equation  [16]  for  the  vibrating 
ellipsoids  described  earlier  have  been  added  to  Figures  3  and  4  where  these 
differ  appreciably  from  the  earlier  results.  The  great  improvement  in  t\e 
accuracy  of  the  results  for  the  higher  modes  is  very  marked,  particularly 
for  the  smaller  length/beam  ratio  ellipsoid.  The  results  with  the  full 
boundary  condition  are  everywhere  almost  identical  to  the  exact  values. 

In  computing  these  values,  the  slopes  b'  of  the  ellipsoids  were 
estimated  numerically  from  the  given  radii  b  at  the  collocation  points  in 
order  to  reduce  the  amount  of  data  needed  to  specify  the  shape  of  each 
ellipsoid  to  that  normal  for  strip-flow  calculations.  Since  this  method 
is  clearly  adequate  for  the  rather  extreme  slopes  involved  in  ellipsoids, 
it  shoul  '  also  be  satisfactory  for  ships  where  shape  changes  are  less 
severe.  Thus  in  applications  to  ship  vibration,  even  the  more  exact  form. 
Equation  [19],  of  the  three-dimensional  flow  approximation  need  involve 
no  more  data  than  presently  necessary  for  the  usual  strip  method. 

The  results  given  by  Equation  [19]  for  ship  vibration  have  also 
been  computed  for  the  destroyer  case  given  earlier.  The  frequencies  com¬ 
puted  are  given  in  the  last  column  of  Table  3  and  show  a  slight  improve¬ 
ment  in  the  predicted  first  mode  frequency.  Again,  the  values  of  b  were 
estimated  numerically  from  the  equivalent  radii  at  the  collocation  points 
so  that  no  extra  data  were  required  in  the  calculation.  As  expected,  the 
improvement  resulting  from  the  use  of  the  exact  boundary  condition. 

Equation  [5],  instead  of  the  approximate  form.  Equation  [6],  was  much  less 
marked  for  the  ship  than  for  the  ellipsoid.  The  mode  shapes  for  the 
destroyer  showed  no  significant  change. 

CONCLUSIONS 

A  method  is  proposed  for  approximating  the  effects  of  the  full 
three-dimensional  flow  around  a  ship  undergoing  transverse  vibration.  It 
provides  an  alternative  to  the  "strip  theory"  usually  used  in  the  calcu- 
lation  of  vibration  frequencies. 
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Although  the  matrixes  involved  are  slightly  more  complicated  to 
set  up,  all  vibration  frequencies  and  mode  shapes  can  be  found  from  the 
eigenvalues  and  vectors  of  a  single  matrix  equation.  In  the  strip  method 
it  is  necessary  to  consider  a  different,  matrix  equation  for  each  mode 
shape.  The  same  data  are  required  for  either  the  three-dimensional 
approximation  or  the  strip  method. 

Where  unusual  mode  shapes  are  involved  (e.g.,  for  damaged  ships 
with  a  very  weak  section)  or  for  the  more  extreme  shapes  of  ship  (small 
length  to  beam  ratios),  the  proposed  method  will  give  better  results  than 
with  strip  theory.  Otherwise,  if  separate  consideration  of  each  mode  is 
acceptable,  strip  theoiy  is  easier  to  apply  and  gives  very  similar 
results. 

The  close  agreement  between  strip  theory  results  and  those  from 
the  three-dimensional  analysis  indicates  that  the  remaining  discrepancies 
between  the  experimental  and  computed  vibration  frequencies  and  shapes  are 
largely  attributable  to  inadequacies  in  the  specification  of  the  elastic 
stiffness  characteristics  of  ships,  such  as  cross-sectional  inertia  and 
shear  area,  rather  than  to  inadequate  representation  of  the  hydrodynamic 
forces . 
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Figure  2  -  Coordinate  System  for 
Vibrating  Ellipsoid 
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FORCE  PER  UNIT  LENGTH 
f(2)/nob2 


(DISTANCE  FROM  CENTER  OF  ELLIPSOID/ (HALF  LENGTH) 

Figure  3  -  Comparison  of  Three-  and  Two-Dimensional  Flow 
Approximations  with  the  Lxact  Solution  for  a  Vibrating 
F.llipsoid  (L/B=iO) 
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(DISTANCE  FROM  CENTER  OF  ELLIPSOI 0/ ( HAL F  LENGTH) 

Figure  4  -  Comparison  of  Three-  and  Two-Dimensional  Flow 
Approximations  with  the  Exact  Solution  for  a  Vibrating 
Ellipsoid  (L/B=5) 
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for  a  Destroyer 


26 


APPENDIX 


EXACT  SOLUTION  FOR  TRANSVERSE  VIBRATION  OF  AN  ELLIPSOID 

OF  REVOLUTION 

The  coordinate  system  used  to  solve  this  problem  is  the  ellipsoidal 
set  (c,  u,  6)  illustrated  in  Figure  2.  The  curves  C  =const  and  u  =const 
represent  confocal  (and  so  orthogonal)  systems  of  ellipsoids  and  hyper¬ 
boloids  of  two  sheets,  respectively.  In  terms  of  the  cylindrical  polar 
coordinates  (r,  6,  z) , 

,  1/2  -  1/2 

r  =  k(4  -1)  Cl-P )  [A-l] 

z  =  k  ;  p 


If  £  =  x,  represents  the  given  ellipsoid,  which  is  assumed  to  be  of  total 

2  2  1/2 

length  2a  and  maximum  diameter  2b,  then  k  =  (a  -  b  )  and 


-1/2 


=  (1  -  8  ) 


where 


8  =  b/a  [A- 2] 

If  the  vertical  velocity  distribution  along  the  axis  of  the  ellipsoid  is 
v(z/a) ,  then  the  boundary  condition  on  the  surface  of  the  ellipsoid  will 
be 
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.  r 

dv  .  .  ,"1 

n 

>n 

■> 

u 

v  cos  X  -  b  sin  A 

cos 

[A-3] 


which  allows  for  the  rotation  of  cross  sections  due  to  bend'.ng.  If  v(z/a) 
is  a  polynomial  of  degree  N,  then  the  solution  for  the  velocity  potential 
can  be  written  (see  Reference  6)  as 


N+l 

0  =  an  P^(y)  cos  9  [A-4] 

n=l 

where  P*(u)  and  Q^U)  are  associated  Legendre  functions  of  the  first  and 
second  kinds,  respectively. 

Substituting  this  expression  for  $  into  Equation  [A-3]  gives 


bk  dv 
a  dy 


N+l 

■z 

n=l 


n  dc 


d  P 
_ r 

dp 
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,  the 


(since  z=au  on  the  ellipse), 
equation  becomes 


Since  also  v(z/a) 


ZV 


n- ' 

Putting  =  -  ^ 
N+l 
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d  Q1 
xn 

ak 

d? 
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C=;o  n. 

d  Q1 
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Z Xn  [V«>  -  Z 

n=l  n=l 


1  -  (n-l) 6  2 

n  vn  u 
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where 


+1 

'an  ’  \  V>d“  *  0 

-1 


n  <  m  or  (n-m) 


_  7  n(n-l)  —  (n-ra+2) _ 

(n+m+1) (n+m-1) . . . (n-m+5)  *  1  '} 


Thus,  given  the  values  v  ,  the  A  and  hence  the  a  are 
’  &  n  n  n 

Kith  the  coefficients  a^  known  in  Equation  [A-4] 
bution  f(z)  on  the  ellipsoid  is  given  by 

7T 

f(z)  =  -  2  J  b(z)  p  cos  9  d6  where  p 
o 


so  that 


N+l 


f(z)  =  -  7rpb(z)  y  an  qJ(?o)  P*(z/a) 


n=\ 


n-l 

u 


to  u  gives 


1,  ....  N  +  1) 


odd 

even 

easily  found. 

,  the  force  distri- 

=  0  6 
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A  short  computer  routine  has  been  written  to  compute  values  of  f(z)/7T0b~, 
given  b/a  and  Vy..,  v5,  for  a  series  of  values  of  z.  For  the  routine,  N 
is  restricted  to  4  since  this  is  sufficient  to  represent  the  vibration 
modes  of  principal  interest. 
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